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Abstract 

We construct the quasi- stationary (QS) distribution for two models of heterogeneous catalysis 
having two absorbing states: the ZGB model for the oxidation of CO, and a version with noninstan- 
taneous reactions. Using a mean-field-like approximation, we study the quasi-stationary surface 
coverages, moment ratios and the lifetime of the QS state. We also derive an improved, consistent 
one-site mean-field theory for the ZGB model. 
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I. INTRODUCTION 

The Ziff-Gullari-Barshard(ZGB) model [l| for the oxidation of carbon monoxide (CO) on 
a catalytic surface exhibits phase transitions between an active steady state and absorbing 
or "poisoned" states, in which the surface is saturated either by oxygen (O) or by CO. This 
model and its variants have stimulated much interest; maiiy studies of the ZGB model have 
been p*...ed us., de.e...Mie .ean-fie.d „ jflQ and Mo,.e Ca.o s„nu.a- 

nnn 

tions. For reviews see p, M,IH-) In the present work we analyze quasi-stationary properties 
of the model. In the process, we revisit the one-site mean-field theory and formulate an 
improved version thereof. 

The surface reactions described by the ZGB and allied models follow the Langmuir- 
Hinshelwood mechanism Q], 

COgas + f COads 
02gas + 2V 20ads 
COads + Oads CO2 + 2V, 

where v denotes a vacant lattice site, and the subscripts indicate the state (gaseous or 
adsorbed) of the molecule. The product CO2 is understood to desorb immediately on for- 
mation. QOgas molecules per site while O2 arrives at rate (1 — Y) = Y. 

Monte Carlo simulations jpy] reveal that the model exhibits two phase transitions. For Y < 
Hi = 0.391 the lattice eventually poisons with O while for Y > 1/2 = 0.5256 it poisons with 
CO. For yi < Y < 1/2 the system attains a reactive steady state. The first transition is found 
to be continuous while the second is strongly discontinuous. Mean-field-like descriptions 
based on cluster approximations reproduce qualitatively the behavior of the system in 
the thermodynamic limit. 

Large but finite systems appear to be stationary on any reasonable time scale (away from 
the phase boundaries), but eventually end up in one of the absorbing states (0-poisoned 
or CO-poisoned), via a fluctuation with a very small, but nonzero probability. Mean- field 
theories ignore such fluctuations, and so are incapable of treating finite systems. Since 
simulations and other numerical methods usually study finite systems, it is useful to develop 
approximate theoretical descriptions that account for finite system size. Experimental study 
of nanoscale catalytic systems is also of current interest |9j . The notion of a quasi-stationary 



(QS) distribution has proved to be a powerful tool in this context. For such models, the QS 
distribution describes the asymptotic (long-time) properties of a finite system conditioned on 



survival 



10 



11 



12l |. The quasi- stationary properties converge to the stationary properties 



when the system size N ^ oo. 

In Ref. jlJ] QS distributions (QSD) are derived for various continuous-time Markov 
processes, e.g., the contact process and the Malthus-Verhulst process. In Ref. |lj| the QSD 
for a discrete-time process is derived, using the Domany-Kinzel cellular automaton as an 
example. In all of these studies, mean-field-like approximations are used to find the QSD. 
Previous studies have focused on continuous phase transitions to an absorbing state, with 
the exception of Schlogl's second model known to exhibit a discontinuous transition 
in mean field approximation, which was shown in Ref. (l^ to exhibit a bimodal QSD at 
the transition. Such a bimodal distribution (e.g., of the order parameter) is a hallmark 
of a discontinuous phase transition. The phase transition in Schlogl's second model, with 
nearest-neighbor interactions, is however known to be continuous in d < 4 dimensions Q|. 
Hence, aside from the particular interest of the ZGB model in the context of catalysis, the 
present work presents a QS analysis of a model known to exhibit a discontinuous transition. 

The balance of this paper is organized as follows. In Sec. II, we briefly review the mean- 
field theory for the ZGB model at the site level, introducing modifications that take into 
account the prohibition against 0-CO nearest-neighbor (NN) pairs on the surface. In Sec. Ill, 
we define the QS probability distribution and show how to derive it in an efficient way. 
This method is used to generate the results reported in Sec. IV; results for the case of 
noninstantaneous reactions (NI model) are also presented. Sec. V contains a summary of 
our results. The Appendix treats the ZGB model on a complete graph. 



II. SITE MEAN-FIELD APPROXIMATIONS 



The simplest mean-field description of the ZGB model is the site approximation. By 
assuming spatial homogeneity and neglecting all correlations between the occupancies of 
different sites, one obtains a closed description in terms of the CO and O coverages, i.e., 
the fraction of sites occupied by C0(0). To simplify the notation, from this point we will 
denote a CO molecule by c, a vacant site by v, and an O atom by o. The coverages will be 
denoted by 9i, for i = c, v or o, with the normalization 9c + 9o + Ov = ^- 



In the ZGB model reactions between adsorbed c and o are instantaneous, so that there 
is a strict prohibition against c-o NN pairs. Previous formulations of the site approximation 







effectively ignored this prohibition, leading to the appearance of an absorbing state 
in which a small fraction (about 0.2%) of the sites are occupied by o, while the rest bear 
c molecules. Elimination of this nonphysical effect is one motivation for reformulating the 
approximation; providing a systematic basis for the QS analysis is another. In the absence 
of any prohibitions, the site approximation treats each site as statistically independent, so 
that the joint probability % for a NN pair of sites to be in states i and j is 6i6j. Our problem 
is to formulate an approximation that treats sites as independent, except for the fact that 
= ^co = 0. (We use 6ij to denote the joint probability that, on the square lattice, the 
leftmost or lower site of the pair is in state i and the other site in state j. By symmetry 

On a square lattice of N sites (with periodic boundaries, so that each site possesses 
exactly four NNs), there are 2N NN pairs. Let Ni be the number of sites in state i and Nij 
the number of NN pairs with states i and j, using the same convention as for 6ij. Then we 
have the relations 

= ^{N,, + N,o + N,,) (1) 
No = ^{Noo + N,o) (2) 



and 



iVc = ^(iVcc + iV,,) (3) 



Letting 6o\o denote the conditional probability that a site bears an o atom, given that a NN 
of this site does, we have 

Noo = 2Noeo\o = ^No-^ (4) 

t>o + Vv 

where we used the fact that if a site bears o, then its NN cannot harbor c. By the same 
reasoning 

No. = 2No-^ (5) 

Oo + 



(6) 



and 



Nr. 



2Nr 



'9c + 9^ 

Combining Eqs. (0), (0) and ((7j) we then have 



2N, 



1 



Or 



(7) 



(8) 



6o + 6u 9c + 

Note that in general iV^ < 2Ny6y, the value it would have in the absence of the prohibition 
against c-o NN pairs. In order to eliminate such pairs, vacancies must, as it were, be 
redistributed to create more v-o and v-c pairs than would be present under completely 
random mixing, and in this redistibution the number of v-v pairs is diminshed. 

Now using 6ij = Nij/2N and 6i\j = 6ij/6j, we may proceed to assign rates to the various 
adsorption/reaction events in the model. Consider, for example, nonreactive adsorption of 
c. This occurs at an intrinsic rate of Y, and requires a vacant site whose four NNs are 
free of o. The probability that a site does not harbor o, given that its NN is vacant, is 
1 — 0o\v = 6v/{6o + Ov), so that the rate Wi (or expected number of events per site and unit 
time) of nonreactive c adsorption is 

Wi = Ye,i^^^^\ =Yd„a (9) 



The rates for the remaining events, listed in Table I, are found in the same manner (for 
further details see . We use the notation 

6o Gc 



Or. 



1 - 



Oo + Oy 9c + 9jj 



and 



7 



(10) 



(11) 



,9c + 9v, 

From the rates listed in Table I, one readily derives the following equations of motion for 
the coverages: 



9r, 



2Y9,,^-Y9,il-a) 



9c = Y9,a-2Y9,,{l--f) 



(12) 



Numerical integration of these equations reveals that for Y < Y^ = 0.56101, starting 
from an empty lattice, the system attains an active steady state. Yg marks a spinodal point. 



the limit of stability of the active phase |7| . The site approximation does not reproduce the 
0-poisoning transition, yielding Y^"^ = 0. 

The present formulation is an improved version of the one-site approach presented in P| . 
In the latter formulation one finds = 0.5615, but, as noted, the final state (for Y > Y^) 
has a small admixture of o in the c-poisoned phase, in contradiction to the prohibition 
against c-o NN pairs. The present formulation eliminates this flaw and provides a slightly 
improved estimate for Y^ . 

In the preceding discussion we consider the thermodynamic limit (there is no reference 
to the system size). The model exhibits an active phase below the spinodal value Y^ . In 
finite systems it is not possible to avoid poisoning. But for large system sizes, we expect, 
for a range of Y, a long-lived metastable state, which we call the quasi-stationary state. 
Increasing the system size, the lifetime of the QS state goes to infinity, and its properties 
become equivalent to that of true active stationary state. In the next section we use the QS 
probability distribution in order to study finite systems. 



III. QUASI-STATIONARY PROBABILITY DISTRIBUTIONS 

The dynamics of a stochastic process such as the ZGB model is governed by a master 
equation (ME) for the probabilities in state space; the latter consists of all possible con- 
figurations of the system. Since the transition rates of the model do not satisfy detailed 
balance, one cannot determine the steady-state probability distribution without solving the 
ME. We may write a simplified ME, based on a cluster approximation, which, depending on 
the truncation level, includes certain fluctuations. 

In the simplest case, one truncates the ME at the site level. This means that the state 
of the system is specified by the site occupation numbers c, o and v, without regard to 
the relative positions of the adsorbed species on the surface. (Evidently, the fluctuations 
preserved at this level are those in the total occupancy.) Consider a square lattice of sites, 
so that c + o + v = N. The state of the ZGB model in the site approximation is given by the 
pair (c, f); the model defines a Markov chain in the c-v plane (see Fig. 1), with 6^ = c/N, 
etc. The 0-poisoned and CO-poisoned (absorbing) states are (0,0) and (A^, 0), respectively. 

In the site approximation, the probability distribution Pc^v{t) follows the ME 



Pc,v{t) = ^ W{^c' ■u')^(c,v)Pc' ,v'{t) — Pc,v{t) ^ W{c,v)^(c' ,v>) (13) 

where W(c,v)^{c',v') is the transition rate. In terms of the ZGB rates shown in Table 1, we 
have: 

5 
i=l 

In Eq. (fT^ it is understood that each of the Wi is to be evaluated for the same values 
as the associated probability factor (that is, in the first term, Wi must be evaluated for 
(c — 1, f + 1), etc.). 

"Poisoning" in the ZGB model corresponds to the system falling into one of the ab- 
sorbing states (see Fig.l). This means that, for a finite system, the stationary probability 
distribution is nonzero only for absorbing configurations, that is, 

lim Pc,v{t) = p5{c,v),(ofi) + (1 - p)S(c,v),iN,o) (15) 

t — ^oo 

(The weight p depends on A^, Y, and the initial distribution.) 

We introduce the QSD via the hypothesis that, as t oo, the probability distribution 
conditioned on survival attains a time-independent form. Then the QS probability distribu- 
tion Qc^y is given by 

Qc, = lim PcAt)/Sit) (16) 

t^oo 

where S{t) is the survival probability, i.e., the probability that the system has not fallen into 
the absorbing state up to time t. In the Appendix, we present results on the QS distribution 
for the ZGB model on a complete graph, which admits an exact description in terms of a 
single variable. 

One method for generating the QSD is by integrating the master equation until the 
distribution P^"*,^ = Pc^y{t)/ S(t), (with S'(t) = 1 — Po,o(i) — -P/v,o(^)) attains a time-independent 
form. For moderately large systems (A^ > 100), standard schemes such as the fourth-order 
Runge-Kutta method require long integration times, because a small time step is needed 
to avoid numerical instabilities. An alternative method is available (2Q]|, based on writing 
the evolution in the form: 



PcAt) = -WcvPcAt) + (17) 

where 

Wc,v = (18) 

and 

^c,v{^) = XI ^ic'y)^ic,v)Pc',v'it) (19) 
i;'>l,c' 

Suppose the probabihty distribution has attained a QS form at some time t. Using 
normahzation, J2v>i cQc,^i^) ~ 1 ^'i- we have 

^c. = (W^c. - ^o)i^c,., (20) 

where 

Ao= J2 Wic'y)^(c,o)Pc,v = -Y. ^-.-(^) (21) 

v'>l,c' '^^IjC 

is the decay rate of the survival probabihty. 
For the ME, Eq. W^,^ = ^J^^ Wi and 

Ac,v = W^lPc-l,.+l + W2Pc,v-l + W^3^c,.+2 + W^4Pc+l,. + W,P,+2,v~2 (22) 

The decay rate is Aq = W^iP/v-i,! + M/3-Po,2- In the quasi-stationary regime we have 

Qc..= „/'^'\ (^>l), (23) 

which suggests the following iteration scheme: 

Pc,^ = ^Pc,. + il-a)-^^, (24) 

where a is a parameter and Ac^y is evaluated using the distribution Pc,v{t) The new distri- 
bution P^^ must be normalized after each iteration. In this way, one can construct the QS 
state from any normalized initial distribution Pc,d(^)- This iterative method exhibits good 
convergence for a = 0.5. 



IV. RESULTS 



We construct the QS distributions for the ZGB model at the one-site level, via both the 
Runge-Kutta method and the iterative scheme, which as expected yield the same QS distri- 
bution. Using a = 0.5, we find that the iterative method converges to the QS distribution 
about 100 times faster than via the RK method using the same precision. This gain in 
efficiency is associated with the small time step required to maintain numerical stability in 
the RK scheme. In Fig. 2, we show a typical QS distribution, Qc,v, for the ZGB model 
near the spinodal point. It is clearly bimodal, with one maximum at c = A^ — 1, near the 
CO-poisoned state, and another, broader one that corresponds the active state. 

Figure 3 shows the time evolution of the coverages (conditioned on survival), for two 
different initial conditions, and Y fixed at the spinodal value obtained from the MF 
theory of Sec. II. It is clear that, although the relaxation to the QS values is very different 
in the two cases, both reach reach the same QS coverages. The relaxations to the QS state 
for Y < Yg^ and Y > Y^ are compared in Fig. 4 (Since the iterative method provides only 
the QS distributions, the evolution is calculated via the Runge-Kutta method in this case.) 

It is also possible to discuss QS properties for Y > Yg, where an active stationary state 
does not exist, even in the thermodynamic limit. This is a finite-size effect, analogous to a 
nonzero magnetization in the Ising model above the critical temperature, on a finite lattice. 
This effect was also reported for the contact process on a complete graph [iS^ . 

The marginal QS distribution (Fig. 5) becomes bimodal as F — > Y^ , consistent with 
a discontinuous phase transition and the emergence of a hysteresis loop when oo. The 

lifetime r of the QS state is given by 

r = 1/Ao. (25) 

As seen in Fig. 6, r becomes very large as we increase A^, for Y < Y^^ . The lifetime increases 
exponentially with A^ (Fig. 6 - inset), as found for the CO poisoning time on finite lattices 



Fig. 7 shows the mean number of adsorbed CO molecules, (c), and of vacant sites (v), in 
the QS state. It also displays the QS moment ratios rric = (C^)/(C)^ and = jV'^) I (V)'^ 
as functions of Y . The quantity m is analogous to Binder's reduced fourth cumulant 



22 



at an equilibrium critical point. At a continuous phase transition to an absorbing state, 



we expect m for the order parameter to attain a universal, size-independent value, due to 
scaling of the probability distribution 



25|. In the present case of a discontinuous phase 



transition, and rriy can be expected to grow due to the bimodal nature of the probability 
distribution. Very large values of m^, can be attained as the coverge 6^ tends to its minimum. 
For the finite systems studied here, the peaks of rric and do not fall at the same 
values of Y, as would be expected for an infinite system. The positions of the maxima 
appear to approach one another as — oo; an extrapolation to TV — >■ oo confirms the peak 
values occur at as given by mean-field theory. 

A. Noninstantaneous reactions 

One of many possible extensions for the ZGB model that have been proposed relaxes the 
condition of instantaneous reactions between adsorbed O and CO. The key difference 
between this NI model and the ZGB model is that there is no longer any prohibition against 
c-o NN pairs. There is now a finite mean reaction time for O and CO adsorbed at NN sites 
to form CO2; the reciprocal of the reaction time is the reaction rate R. (We continue to 
assume instantaneous desorption of CO2.) The transitions are listed in Table II; rates are 
obtained as in Table I. (The prefactor 4 in rate W3 is the coordination number of the square 
lattice.) 

In this case, the mean-field equations are: 

60 = 2(1 - Y)e,^ - AROoe, 

= Ye, - AROJ, (26) 

A simple analysis of equation ()26|) shows that Y^ ~ \/2R as i? — 0, grows monotonically 
with i?, attaining the limiting value Y^ = 2/3 as i? ^ 00. This value is different from the 
spinodal found in the site-approximation for the ZGB model, which involves instantaneous 
reactions. Since lim/j_,oo 7^ ^^zgbi limit of infinite reaction rate is singular. The 
singular nature of this limit derives from the fact that the instantaneous reaction condition 
imposes certain constraints on configurations, thereby altering the algebraic form of the 

n 

mean-field equations (a more extensive discussion may be found in ^j). In particular, for 
R 00, 9o and 9c cannot both be nonzero, whereas both coverages are nonzero in the active 
phase of the ZGB model. 



The probability distribution for the NI model follows the master equation: 

Pc,y{t) = RlPc-l,v+l + R2Pc,v+l + Rs [Pc+l,v-l + Pc,v-l] - [Rl + R2 + 2i?3] PcAi)- (27) 

We turn now to a discussion of the QS properties, obtained using the iterative method 
based on Eq. (j23p . The lifetime of the QS state, in the reactive window, increases with R 
(see Fig. 8). The displacement of the maximum position of r (as a function of Y) parallels 
the variation in Y/'{R). The maximum lifetime (i.e., varying Y) tends to a finite value, 
different from that from the original ZGB model, when R 00. This is because the flux 
to the absorbing configurations, which determines the lifetime of the QS state, and which, 
in the site approximation, is given by Aq = WiQn~i,i + W^3Qo,2, is reduced as R increases, 
since the latter reduces the QS probabilities Qjy^i^i and (5o,2- 

Figs. 9 and 10 show the QS coverages 6^ and 6^. As R increases the transition becomes 
progressively "harder," i.e., the coverage 6c just below the transition approaches zero while 
that just above approaches unity. In the limit R —>■ 00, the c coverage is a step-function: 
9c = Q{Y — Y^). As before, the QS coverages tend to the macroscopic curves (furnished by 
mean- field theory) when N ^ 00. Due to the low c coverage in the active phase for large 
R, the moment ratios differ strongly from those of the ZGB model. But as always, the 
peak in nic tends to the spinodal as iV — > cx) (see Fig. 11). 



V. SUMMARY 



We study the quasi-stationary properties of the ZGB model for the catalytic oxydation 
of CO, and of a related model with noninstantaneous reactions. The QS distribution is 

n 

constructed numerically via a bivariate extension of the iterative method [20|. For the 
simple case of a complete graph, we obtain exact (numerical) results for QS properties. 
Although our analysis of the models with NN reactions on a square lattice is based on 
site-approximations, which cannot capture critical fluctuations, the approach allows one to 
put some fluctuations "back into" mean-field theory, specifically, fluctuations in the overall 
coverages, since we assume spatial homogeneity. This enables us to investigate finite-size 
effects and moment ratios that are out of reach with other approaches. We derive a revised 
site approximation for the ZGB model that respects the prohibition against c-o NN pairs, 
eliminating the inconsistency noted in previous formulations 0, Q|, while providing a 



reasonably good estimate of the spinodal point. 

The present study confirms two fundamental properties of the ZGB and NI models: (1) 
the discontinuous phase transition is signaled by a bimodal QS probability distribution; (2) 
the lifetime of the QS state grows exponentially with system size, in the reactive phase. The 
latter implies that the system will typically spend a long time in the QS state, following 
a relatively brief transient, except near the spinodal, where relaxation slows down. We 
have also noted sharply increased values of the moment ratios rric and m„ in the vicinity 
of the transition. Since the CO2 production rate is closely related to 9^, one should expect 
to observe large fiuctuations in product concentrations in the gas phase, as the spinodal 
is approached. One may envision using these fiuctuations as part of a control scheme that 
adjusts the CO partial pressure (represented by the parameter Y in the models), to optimize 
production while avoiding poisoning. 

An interesting result is that at the level of mean-field theory, the ZGB model is not 
equivalent to the infinite-reaction-rate limit of the model with noninstantaneous reactions. 
The ZGB model is characterized by the constraint that any c-o NN pair react before any 
further adsorption events, a condition that is also realized in the NI model, for sufficiently 
large R. Mean-field theory is however incapable of representing this situation, since, in 
the absence of a prohibition against c-o pairs, the NI model in the R —* 00 limit simply 
prohibits simultaneous nonzero coverages of O and CO. Since truly instantaneous reactions 
do not exist, one might suppose that the NI model, with an appropriately large rate R, 
provides the correct description of fast (but not instantaneous) surface reactions. In mean- 
field approximation however, the large- i? limit of the NI model appears less realistic than the 
ZGB model, which permits simultaneous presence of adsorbed O and CO. In this sense, the 
choice of model in the description of cooperative kinetics should involve the approximation 
that one intends to use, not only the intrinsic nature of the model. 

QS analysis promisses to be a useful technique in the study of catalysis models on nanome- 
ter scales pj. Possible directions for future work include the study of metastability via QS 
distributions, and Monte Carlo simulations of the QS state, in catalytic models and other 
processes exhibiting absorbing states. 

ACKNOWLEDGMENT 

This work was supported by CNPq and CAPES, Brazil. 



APPENDIX: ZGB MODEL ON A COMPLETE GRAPH 



A complete graph is one in which all sites are neighbors. When formulated on such 
a structure, a stochastic lattice model defines a stochastic process defined by one or a 
'ew variables, thereby permitting an exact analysis, as in the case of the contact process 
isl ]. Since each site interacts with all others, the behavior is mean-field-like. Due to the 
prohibition against c — o NN pairs, in the ZGB model on a complete graph, only one species 
(c or o) may be present at any moment. Letting n = c — o, we can describe a system of 
sites completely by a single variable, with n = —N representing the 0-poisoned state and 
n = N the CO-poisoned state. The number of vacant sites is f = — |n|. 

Each adsorption event leads either to n^n + loT to n^n — 2. The transition rates 

are 

Wn-.n+l = y(i-^-^) (28) 



and 

iy„ n = y M - ... 

N J \ N 

The process obeys the master equation: 



P„(t) = YUn-lPn~lit) + YUn+2 { Un+2 " ^ ) Pn+2{t) 



-Ur. 



Y + Y\un-^ 



N ^ 

Pnit), (30) 



where m„ = ^1 — The macroscopic equation for this process 21 1 is given by: 

0(t)=F(l-|0(t)|)-2F(l-|0(t)|)l (31) 

where limjv^oo0(^) = n{t)/N. In this limit, the system poisons with CO when Y > 2/3. 
For Y < 2/3 there is a nontrivial fixed point at 0* = 1 — (1 — Y)/2Y witch is the attractor 
if 0(0) > —(f)*. The point 0* is the reactive stationary state for an infinite system. 

Analyzing the master equation via the iterative method described in Sec. Ill, we obtain 
the QSD shown in Fig. 12. The QSD changes qualitatively as Y approaches the spinodal. 
For small Y, the distribution has a single maximum near 0*. For y ^ 2/3 the distribution 
becomes bimodal, with one peak near 0*, related to the active phase, and another at n = A^— 
1, related to the CO-poisoned state. Fig. 13 shows the QS activity for this model; as oo, 
it tends to the macroscopic limit (dashed line). As verified for the site approximation, the 



lifetime grows exponentially with N in the active region (0 < F < 2/3) and linearly in the 
nonreactive region. 
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TABLE I: Site mean-field theory for the ZGB model 



Event 


Rate 




{c^v) to 


V + c |— > c 


Wi = 




{c + l,v-l) 


V + c [-^ V 


W2 = 


Ye,{l-a) 


{c,v + 1) 


vv + oo j— oo 


w^ = 


YOy^^^ 


{c,v-2) 


vv + oo j— ^ vo 


Wa = 


2F^^„7(1 - 7) 


{c-\,v) 


vv + oo j— > vv 


W5 = 


F^^„(l-7)2 


(c-2,v + 2) 



TABLE IL Site mean-field theory for the finite reaction-rate model 
Event Rate (c, v) to 

ci Wi = e^Y {c+l,v-l) 
00 i W2 = 2elY {c,v-2) 
oc T W3 = mOcOo {c-\,v + 2) 



FIGURE CAPTIONS 



FIG. 1. Transitions in the ZGB model (site approximation). Circles denote absorbing states. 

FIG. 2. A typical bimodal QS distribution Qs{c,v) near the spinodal point (Y — 0.56101). 

FIG. 3. Time evolution of the mean values conditioned on survival for Y — 0.56101 and 
system size N — 64. Initial conditions, from top to bottom: half c and half v; full v; half a 
and half v. Dashed lines show the respective unconditioned values. 

FIG.4.Relaxation for the QS mean values of c for y = 0.50, Y = 0.56101 and Y = 0.60, 
from bottom to top. A" = 64. Dashed lines: unconditioned mean values 

FIG. 5. Marginal distribution Qc for Y=.54, .55, .56, .57, .58. 

FIG. 6. Lifetime r of the QS state as function of CO arrival rate Y, for the ZGB model. 
System sizes A = 36, 64, 100, 256, 400, 512 from bottom to top. Inset: Lifetime as function 
of the system size A for y = 0.50. 

FIG. 7. Mean values conditioned on survival (c), (v), and respective moment ra- 
tios, rric and rUy for the original ZGB model without desorption. System size A = 
16,36,64,100,256,400,700. 

FIG. 8. Lifetime r of the QS state for the NI model. System sizes A = 100 (lower sohd hues); 
A = 256 (dashed hnes); A = 400 (upper sohd lines). Reaction rates R = 0.1, 1, 10, andlO^ 
from left to right. The dashed line is the lifetime for the ZGB model. Inset: detail of the 
same graph with A — 100, showing the lifetime for the ZGB model (dashed line). 

FIG. 9. QS coverage 6c for the NI model; reaction rates as in Fig. 8. System sizes A = 100 
(solid lines) and A = 400 (dashed lines). 

FIG. 10. QS vacancy fraction 6^ for the NI model, as in Fig. 9. 

FIG. 11. Moment ratio rUc for NI model. Reaction rates and system sizes as per Fig. 9. 
FIG. 12. QSD for the ZGB model on a complete graph of 100 sites. 



FIG. 13. Quasistationary activity for the ZGB on a complete graph. The dashed line 
represents the macroscopic hmit, Eq. ()31|) . 
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